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ABSTRACT 

Context. We search for new variable B-type pulsators in the CoRoT data assembled primarily for planet detection, as part of CoRoT's 
Additional Programme. 

Aims. We aim to explore the properties of newly discovered B-type pulsators from the uninterrupted CoRoT space-based photometry 
and to compare them with known members of theySCep and slowly pulsating B star (SPB) classes. 

Methods. We developed automated data analysis tools that include algorithms for jump correction, light-curve detrending, frequency 
detection, frequency combination search, and for frequency and period spacing searches. 

Results. Besides numerous new, classical, slowly pulsating B stars, we find evidence for a new class of low-amplitude B-type pul- 
sators between the SPB and 6 Set instability strips, with a very broad range of frequencies and low amplitudes, as well as several 
slowly pulsating B stars with residual excess power at frequencies typically a factor three above their expected g-mode frequencies. 
Conclusions. The frequency data we obtained for numerous new B-type pulsators represent an appropriate starting point for fur- 
ther theoretical analyses of these stars, once their effective temperature, gravity, rotation velocity, and abundances will be derived 
spectroscopically in the framework of an ongoing FLAMES survey at the VLT. 

Key words. Methods: data analysis; Stars: oscillations; Stars: variables 
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, 1. Introduction The study of stellar oscillations from the ex oplanet data of th e 

mission is part of the Additional Programme d Weiss et al.l2004l) . 
In the forthcoming months and years, the French-European In this paper, we focus on CoRoT's Initial Run in the anticen- 
Space mission CoRoT (Convection, Rotation and planetary tre dire ction of the Galaxy (IRaOl). We first describe our meth- 
Transits |Auvergne et al. || 2009|) , will measure brightness varia- ods of analysiSi as future runs wi n a i so be addressed with the 
tions of tens of thousands of different stars in its exoplanet pro- same met hodology. During IRaOl, which lasted some 55 days, 
gramme, with a precision typically two orders of magnitude bet- the CoRoX satellite was pointed at a field in the galactic planei 
ter than we can achieve from ground-based observations today. towards the galactic anticentre (a » 06 h 4 4 m, 6 » -01° 12'). 
This means that we are able to extend the sample of known pul- We ^ focusing on the stars (e g . Lstankov & Han dler 

sating stars to class members with fainter apparent magnitude. [2005!) and slowly pulsating B stars (SPBs, Waelkensl ll99lT) 

discovered in this field. A short descriptive overview of the 



main properties of these two classes is listed in TableQ] There 
is a small overlap between their instability regions, where hy- 
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Table 1. Overview of the main properties of (3 Cep and SPB pul- 
sators. 
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Fig. 2. Averaged Fourier periodogram of 100 non-variable stars. 
The inset is a zoom on the low-frequency region. A low- 
frequency powerlaw trend is present, as are the 1 and 2 day alias 
and the orbital frequency with many of its harmonics and day 
aliases. 



opacity bump at T ~ 2 x 10 5 K, which means that the excitation 
of modes in those stars is sensitive to the behaviour of opacity in 
the stellar interior and to the abundance of Fe-group elements. 

It would be a huge loss of resources to start analysing all 
stars in the CoRoT database only to study the B-type pul- 
sators. Therefore, we resort to the results obtained in the frame- 
work of the variability classification programme described in 
Debossch er et al.l (|2009), hereafter termed CVC, which extracts 
a basic number of relevant parameters from each light curve to 
compare them with training sets of stars of known type. 

In this paper, we focus on all light curves given a label 
of either £1 Cep or an SPB star. We also include the 6 Set and 
y Doradus stars to limit the red edge of the SPB instability strip 
and to discuss possible overlap between /3 Cep and 6 Set stars, 
on the one hand, and between SPB and y Doradus stars on the 
other, in terms of light curve properties (e.g., frequency values). 
In total, 540 light curves are extracted. We emphasise that the 
results from the CVC classifier are only used to select the initial 
sample of stars from the CoRoT database. 

2. Description of the dataset and pre-processing 

All datasets are taken over a period of 55 d. Two time sam- 
plings occur, one measurement every ~ 32 s or one every ~ 5 12 s. 
This sampling is not always constant for one light curve, it may 
change from low temporal resolution to high temporal resolu- 
tion during the observing run. The magnitudes of the stars in the 
sample lie between 12.1 and 16.6 and vary over different scales 
(Fig.0. 

To have an idea of the different kinds of instrumental effects 
we can expect, as well as their impact on the analysis, we anal- 
ysed 100 constant stars and averaged their Fourier periodogram. 
They were selected with the criterion that their dominant fre- 
quency has a p-value >0.01. From Fig. [2] we can identify at 
least three different groups of features: a low-frequency power- 
law trend due to instrumental drifts (red noise), which becomes 
negligible around0.4d _1 , peaks at 1, 2, and4d _1 and peaks con- 
nected with the orbital frequency due to temperature changes. 
We refer to lAuvergne et al.l (120091) for a discussion of the ori- 
gin of these instrumental peaks and to Sa madi et al.l (120061) for 



a description of the CoRoT data pipeline and treatment of data 
during the passage of the South Atlantic Anomaly. 

The most variable feature is the low-frequency powerlaw 
trend, because of which we introduced a pre-processing step 
to minimise its potential influence. Two phenomena add up to 
this trend: discontinuities in the light curve and long-term ef- 
fects, possibly of instrumental origin. It is impossible to decide 
whether or not these longterm trends are real variations in the 
star's light curve, but because we do encounter trends suspi- 
ciously linear or exponential, we choose to eliminate them. This 
also means that we effectively remove possible oscillating be- 
haviour and intrinsic trends on the timescale of the time series 
itself. 

Long-term trends are removed in the following iterative way. 
First, jumps are detected by sliding two adjacent bins over the 
lightcurve. Physically, we only expect to find a gradual, con- 
tinuous rise or decay in the measured flux. To accomodate for 
the dependence of amplitude and present pulsations, we com- 
pare the difference between the average of the bins to the overall 
difference. If at any point the difference is too high to be consid- 
ered continuous, as defined by the entire light curve, we identify 
the point as a discontinuity, perform a piecewise detrending of 
the two separate parts of the lightcurve and try to detect other 
jumps. If no more jumps are detected, the piecewise detrending 
is redone using the original light curve and the discontinuities 
previously detected. There is some arbitrariness in the choice 
of parameters (bin width, threshold value etc.), so they were 
fixed empirically. For the piecewise detrending, the best fit of 
three simple models was chosen for each part; either a linear, 
a quadratic, or an exponential trend. An example is provided in 
Fig.S 

3. Frequency analysis methodology 

3.1. Basic treatment 

To analyse all jump-corrected and detrended light curves, we de- 
scribe the flux F(f) as 

F(ti) = ii + C + J] A, sm(2n(fit + 0,)) + S. (1) 
i=i 
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Fig. 3. Correction of one discontinuity in the light curve of SPB 
candidate 102815542. The first part and second parts were de- 
trended using a quadratic polynomial and exponential decay 
model, respectively. 



Fig. 4. Result of a simulation study to determine the number of 
independent frequencies in the CoRoT datasets {sold line: evenly 
spaced case, Nj = 2N b s , dotted line: fit through the simulation 
results). 



where e, is a Gaussian noise component, /i the observed av- 
erage flux level, and C a usually small correction parameter 
on yi. In this model, there are n nonlinear parameters for 
which the values are estimated usi ng the Lomb - Scargle versio n 
of the discrete Fourier transform dLombl[T976t IScargldll982l) . 
The linear parameters A, and 0, describe the amplitude and the 
phase of the corresponding pulsation, respectively, and are es- 
timated by ordinary least squares regression. The error on the 
frequency determination (and other parameters) is calculated ac- 
cording to Montgo mery & Q'Donoghueldl999h . These error de- 
terminations underestimate the real error because of correlation 
effects in the data. In our sample, this factor lies between 1 and 4 
for the first frequency, and dec reases for subsequent frequencies 
dSchwarzen berg-Cze rny|200"3l) . 

Traditionally, the frequency acceptance criterion of 
Breger et al] d 19931) . using a signal-to-noise ratio (SNR) lower 
limit of SNR — 4 in the amplitude periodogram for a peak to 
be considered genuine, has been proven to be very successful, 
since very few frequencies identified this way (if at all) had 
to be refuted afterwards. However, there is no rule about how 
large one has to choose the interval in the periodogram to 
use for noise calculation. Second, this method is not suited to 
detecting regions of power excess, either due to densely packed 
frequencies or stochastically excited frequencies. Third, related 
to the previous argument, the SNR method does not guarantee 
that, if a peak is considered as due to noise, there is no other, 
possibly even smaller, peak in the periodogram that is not noise. 
Therefore, the SNR method is not fully objective in its use. 

Instead of using the classical SNR criterion, we make an as- 
sumption on what the datasets should look like when devoid of 
instrumental effects and variability from the star itself. A rea- 
sonable (yet not perfect) assumption is that we are dealing with 
white Gaussian noise and that the number of datapoints is large. 
Degr oote et al] (120091) provide an extensive discussion of the de- 
viation from white Gaussian noise and have shown that it can be 
ignored in the frequency analysis, as long as a correction factor 
for correlated data is used. Under this assumption, the distribu- 
tion of the Scargle power spectrum (normalised with the total 
varian ce of the data) belongs to the exponential family dScargld 
119821) . 

P(z) = exp(-z), (2) 



and thus the noise level is located at a height of z — 1 . In reality, 
we never test a single frequency, but cover a wide range of possi- 
ble frequencies, between and the highest detectable frequency, 
the Nyquist frequency f Ny . In t heory, this number is we ll-defined 
regardless of the time spacing (lEver & B artholdi 1999). In prac- 
tice, however, thi s frequency is usually too high to be of any use 
(e.g. |Koen 2006). A better practical guess for the Nyquist fre- 
quency in uneven time sampling is obtained via 

fyy = ^ O) 

with At in our case the value of the larger timesteps. 

To account for the simultaneous testing of all frequenc ies 
up to fyy, we apply the Bonferroni correction dScargldll982h to 
Eq.©: 

/ , B (z) = l-[l-exp( Z )f'. (4) 

For equally spaced data, theory predicts that the number of inde- 
pendent frequencies iVj between en fy y scales with the number 
of observations N b s as 

Ni = 2N obs . (5) 

In the case of un evenly spaced data, this parameters Ni loses 
its meaning (e.g. Fres cura et al]|2008l) . However, a simulation 
study shows that expression (0, combined with definition ©, 
turns out to be a useful and conservative approximation of the 
true number of independent frequencies. For eight different time 
templates taken from a CoRoT light curve, we generated 10000 
light curves consisting of white Gaussian noise. Next, we cal- 
culated the height z of the maximum peak in the Scargle peri- 
odogram. These were combined, and a false alarm function of 
the form (0| fitted through the data, to determine the free pa- 
rameter Ni. The result can be seen in Fig. [4] and shows that the 
CoRoT data allow the use of Eq. (0. 

Although the sampling rate is usually much higher than 
strictly needed, we choose not to rebin the data or to convert the 
flux to magnitude, because such manipulations alter the shape of 
a light curve. 
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3.2. Search for combination frequencies 

In particular cases, model (Q~|i can be expanded further to lower 
the number of dependent parameters and thus simplify the model 
in terms of fi. Neglecting this effect implies imposing a linear 
model on all observed light curves, thus treating every found 
pulsation frequency as an independent mode. Identifying depen- 
dent frequencies can be seen as a basic correction of the linear 
model for nonlinear effects: 



0.15 



F(t) =n + C +Yj A k sinlXAf/ + <Pk)] 

k=\ 

+ Y^Aisin[2n(fiti + (f>,)] 



(6) 



with 



where iif < n and // represent combination frequencies. 
This might come from nonlinear effects, e.g., nonlinear cou- 
pling between modes dBuchler et alJll997h Handler et alj|2006t 
Degr oote et alj|2009l e.g.). In model ©, the linear parameters 
still have to be fitted, but are sometimes expected to be in phase 
or in anti-phase with each other. In contrast to the window fre- 
quencies, combination frequencies and harmonics are not artifi- 
cially introduced, because they have a real physical interpreta- 
tion. It is thus important to derive them. 

The search for combination frequencies is done in the fol- 
lowing way: for each frequency, the harmonics are identified us- 
ing all empirically found frequencies with higher amplitude than 
the one under investigation. Those are then excluded from the 
list of independent frequencies. For the remaining frequencies, 
all second and third order combinations are tested, again using 
only the frequencies with higher amplitude, and not part of a 
combination so far (in order to avoid higher order 'combinations 
of combinations'). To keep the parameter space manageable, fre- 
quencies of order n are only considered if at least one of order 
in— 1) is previously identified. We are limited by the total time of 
observations, so there is some ambiguity in the choice of possi- 
ble deviation from the exact combination. Because the Rayleigh 
limit of Lr — l/T is a natural measure of maximum uncertainty 
on the frequency determination, we consider this a reasonable 
value to test the presence of combinations. 



3.3. Search for period and frequency spacings 

The variability detected in SPB stars is interpreted as stemming 
from low degree, high-order modes whose periods are known to 
be sensitive to the structure of the stell ar core. The fir st-order 
asymptotic approximation developed by TassouJ (Il980h shows 
that the period spacing (AP) between modes of consecutive order 
and same degree is constant and can be approximated as 



AP = 2jt f L f —dx) 

Ja x 



(7) 



where L = \l(t + l)] 1 ^ 2 (with i the mode degree), x the normal- 
ized radius, xq corresponds to the boundary of the convective 
core and N the Brun t-Vaisala frequency. Moreover, as shown by 
Mi gr ip et~aT] d2008l) . deviations from a uniform period spacing 
are very sensitive probes of the chemical composition gradient 
in the star that develops near the edge of the convective core. A 
sound detection of period spacing in SPB stars would therefore 




4.0 M„ 



Z = 0.02 

i ... i i 



4.3 4.2 4,1 4 39 

Fig. 5. Theoretically predicted asymptotic i — 1 period spacing 
of high-order g modes (AP) as a function of the effective temper- 
ature. AP is computed for main-sequence models with masses 
between 2.5 and 8 M@, with (full line) and without (dashed line) 
overshooting from the convective core. The heavy-element mass 
fraction assumed in the models is Z = 0.02. 



allow direct inferences on the near-core structure of B stars. A 
theoretical estimate of the expected average AP for t = 1 modes 
is reported in Fig.|5j where main-sequence models between 2.5 
and 8 M Q , computed with and without overshooting, are con- 
sidered. We have to recall, however, that such a regular period 
spacing could be strongly affected by the effects of rotation on 
such long oscillation periods (see e.g. Dziembowski et alJI 19931 
for a discussion). 

Analogously to the case of g-mode periods, a regular fre- 
quency spacing can occur for p modes. Other situations where a 
form of spacing can be found, is in frequency multiplets origi- 
nating from stellar rotation, where deviations from the equidis- 
tance case can also occur. A systematic search for recurrent fre- 
quency/period spacings therefore represents a valuable tool for a 
theoretical interpretation of the large number and variety of pul- 
sating stars detected by CoRoT. The total time span of 55 days 
for the Initial Run puts an upper limit of ~ 0.02 d _1 on the fre- 
quency resolution. For rich SPB pulsators, the frequencies are 
expected to be densely packed, thus detecting period spacings 
may be difficult. 

To detect spacings in a list of frequencies ( f\, . . . , /„) or 
periods (Pi, . . . , P„), we use the Kolmogorov-S mirnov test of 
Kawaleri (l988) and the inverse variance test of lO'Dorio ghue 
(119941) . In addition, a third new test is proposed here. For the 
latter, we first construct a vector containing all possible differ- 
ences between the periods. This way, we are left with a vector 
of length Binom(«,2), and elements Ay = |P, - Pj\. For a test 
spacing AP, we calculate for each element of these entries A (/ , 



AP 



AP, 



(8) 



where |_*J is the greatest integer smaller than x. We count the 
number s of elements in the intervals 

s -> [0,e] U [AP-e,AP]- 
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This number s captures the observed extreme values, which are 
defined through the customizable tolerance parameter e. If we 
assume that these differences are uniformly distributed, we can 
express the expected number s e of extreme values by 



2e 
AP 



Hi. 



Application of the plain matrix test implies computation of the 
p value according to 

P(X,e[0,e]U[AP-e,AP]) = ^. 

This test is appropriate if some periods are missing or not de- 
tected. However, we may prefer to only take those pairs into 
account for which the difference between the two members is 
exactly AP, and not a multiple. This we can do using the distinct 
number test. It works the same as the matrix test, except that we 
keep track of the term 



AP 



AP 



from Eq. (8). Only those spacings where m = 1 and r < e or 
where m = and AP - r < e are considered as valid period- 
spacing candidates. The possible disadvantage of using this test 
is that we can end up with a lot of periods belonging to the same 
AP, but where no single triplet is detected. Of course, this can 
be an advantage as well, e.g., in the search for rotational split- 
tings; the matrix distinct number test will detect doublets with 
the same spacing, far apart around different frequencies. Finally, 
if we are interested in the number of periods belonging to each 
test spacing AP, we can simply count the number of different 
periods belonging to each test spacing period set. 

To illustrate the performance of these tests, we simulate a 
collection of periods, in which we have put 3 different spacings 
(in arbitrary units): AP[ = 0.1315 (8 periods), AP 2 = 0.4332 
(6 periods), and AP 3 = 0.3579 (3 periods). On top of that, 
we added 8 random periods, to arrive at a total number of 25 
periods. We additionally allowed each period to deviate from 
the rigid spacing according to a random normal fluctuation of 
cr = 0.0005. The results of the different tests are shown in 
Fig. [6] the Kolmogorov-Smirnov test locates only the most obvi- 
ous period spacing, but suffer s from a strong bifurcation artifact, 
already noticed by iKawalerl (Il988l) . The inverse variance test 
shows lack of power in this dense pack of periods; it is only able 
to detect a spacing of APi/2. This contrasts to the matrix tests; 
their statistic is almost devoid of noise, and only show peaks at 
the true period spacings and their multiples. The matrix test and 
distinct number test show great similarity in their peak structure, 
because they have the same basic calculation method. The rea- 
son the matrix test favours the larger spacing AP2, although it 
has fewer members than APi, is how we set up the statistic: we 
gave larger spacings less chance of happening. 

Regardless of the succes of these tests, we must be careful 
in interpreting the results: because of the occurrence of other 
random periods, we may by chance pick up other pairs of periods 
with the same spacing, but totally unrelated. 

4. Evaluation of the new pulsators 

4. 1 . Additional Stromgren photometry 

To determine the effective temperature T e g and gravity logg 
of the stars in the initial run, Stromgren photometry was ob- 
tained using the Wide Field Camera attached to the Isaac Newton 
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Fig. 6. Outcome of different tests for period spacings on a sim- 
ulated set of periods. Input spacings AP, and detected spac- 
ings APJ were equal for APi = 0.1315 and AP 2 = 0.4332. 
For the third spacing AP3 = 0.3579, the detected spacing was 
AP3 = 0.3576. The (top panel) Kolmogorov-Smirnov test with 
the typical bifurcation, and inverse variance test (second panel) 
do not detect all spacings. The distinct number test (third panel) 
and matrix test (bottom panel) are more suitable to detecting the 
spacings. 



Telescope at the La Palma observatory, Spain. Because of time 
constraints, no (3 index measurements could be obtained. This 
implied that we did not have a direct measurement of the gravity. 
We thus proceeded in an iterative way. First, we computed the 
dereddened photometric indices by using an appropriate spec- 
tral type for each of the pulsators: B2V for the candidate ySCep 
stars, B5V for the candidate SPBs, A2V for the 5 Set candidates, 
and F2V for the y Dor candidates. Next, we used cq to predict the 
ft index following the r elation between these two quantiti es de- 
rived b y lBalonal d 19941) for main sequence stars. In fact. iBalonal 
(1994) has shown that this procedure has advantages over us- 
ing the measured indices to compute ft, due to imperfections in 
the calibrations. With these dereddened indices and ft estimates, 
the effective temperatures and gravities were subsequently esti- 
mated, using Balona's (1994) calibration. 
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The procedure described above delivers an estimate of T e g 
and logg with an internal error, i.e., an error assuming that the 
dereddened input colours are error-free. It is well known, how- 
ever, that this procedure for determining the fundamental param- 
eters suffers from systematic uncertainties connected with the 
limitations of the model grids used to explain the behaviour of 
standard stars. To get a han dle on these systematic err ors, we 
also used the calibrations of iMoon & Dworetskvl d 1985b . which 
are valid along the entire main sequence as well. The standard 
deviation between the values for T e s and \ogg derived from 
Balona's (1994) and Moon & Dworetsky's (1985) calibrations 
was considered to be a good approximation of the systematic 
uncertainties, so these were added to the statistical errors. 

As a test of the validity of the obtained fundamental parame- 
ters for the new B-type candidate pulsators, we placed all stars in 
the HRD and compared their position with the instability strips 
(Fig-0- The numerous new 6 Set stars match the classical insta- 
bility strip perfectly, which gives us confidence that the estimates 
of T e s and log g obtained from our procedure are appropriate. 

From Fig. [7] we see that the classification is very power- 
ful for extracting candidate new pulsators from a large sam- 
ple on the basis of wh ite-light photometry alone (as done in 
Debossch er et al.l d2009l) ). but it is not perfectly reliable on the 
scale of individual stars. Moreover, we do not expect so many 
candidate ySCep stars compared to SPBs and 8 Set stars, consid- 
ering the mass-dependence of the initial mass function and of 
main-sequence lifetimes: this implies that several of them might 
have a different character. We therefore decided to regroup all 
new candidate B stars in the sample, by using the and log g 
information, along with the obtained frequency spectra. This ad- 
ditional classification should allow a better distinction between 
the classes of pulsators along the main sequence than just using 
the white-light photometry alone, since we add information on 
the fundamental parameters. The new clusters we obtain in this 
way are discussed by means of typical and atypical examples, 
some of which are shown in the text, while others are reported 
in the electronic Appendix A. 

4.2. Candidate pulsators situated between the SPB- and 
6 Set instability strips 

The most obvious result from Fig.|7]is undoubtedly the appear- 
ance of p-mode type variations (in terms of frequency values), 
not only far from the theoretical /?Cep instability strip, but also 
outside any instability strip. Part of this discrepancy can be at- 
tributed to the CVC not containing classes of spotted or dif- 
ferentially rotating stars, or ellipsoidal variables. Unavoidably, 
the main frequency overlaps with the /3 Cep pulsation range for 
some of these stars. Moreover, many hot stars in our sample turn 
out to be good Be star candidates (Neiner, private communi- 
cation). These are or will be studied and p ublished elsewhere 
dGufierrez-Soto et al.l 120081; lEmilio et all I2009L and future pa- 
pers). 

The only remaining /?Cep candidate is CoRoT 102813271, 
where two frequencies are detected in the typical p-mode fre- 
quency regime: f\ a 5.79d~', fa « 5.51 d~' and a third fre- 
quency fa « 0.23 d~' in the g-mode regime. Because of its low 
amplitude and long period, the g-mode frequency can possibly 
be influenced by instrumental effects. 

The objects on the red side of the instability strips do 
not constitute one uniform group. That they all vary over the 
same timescale implies that the CVC classified them as /JCep 
stars. However, looking beyond the first few frequencies reveals 
their distinct characteristics: CoRoT 102833548, 102848985, 




4.5 4.4 4.3 4.2 4.1 4.0 3.9 3.8 3.7 
lo gio r eff 

Fig. 7. Location of targets in a (r e ff,logg) diagram. The SPB 
star candidates (white dots, 66 stars) occupy the theoretical in- 
stability strips well (grey dashed lines), except for a few out- 
liers (see text). The same holds for the 92 6 Set candidates (light 
grey dots), but not for the 43 ySCep candidates (dark grey dots). 
Arrows denote outliers: the cool outlier CoRoT 102818325 from 
the ySCep candidates is probably a binary. The two cool out- 
liers on the 8 Set side each have only two detected frequencies: 
CoRoT 102921009 has f x « 8d ',/ 2 = fi/2 (possible a binary), 
and CoRoT 102787451 has fi * 17.4d _1 ,/ 2 ~ 16.5 d ! . 



102922479, and to a lesser extent also the almost constant 
star CoRoT 102850576, have a well-structured frequency spec- 
trum consisting of equally spaced frequencies (Fig. IA.ll in the 
Appendix). The results of the spacing tests and a time-frequency 
analysis of the detrended light curves are shown in Figs[8] and 
[9] and indicate that a spacing of ~ 0.70 d~' is common. The 
amplitude of the main peaks are of the order of 0.1%. The fre- 
quency range overlaps that of /3 Cep stars on the low side, but 
also reaches lower frequencies. The determined temperatures 
place all these stars consistently between the SPB and S Set in- 
stability strips. 

CoRoT 102861067, 102790063, 102790331, 102729531, 
102862454 and 102771057 constitute a second well-separated 
group. They share a complex frequency spectrum containing 
between ~20 and ~170 significant peaks spread over a wide 
range and amplitudes below 0.1%. Two examples are shown in 
Figs |A.2| and |A.3l while schematic oscillation spectra are shown 
in Fig. IA.4l in the Appendix. Dense 'forests' of closely spaced 
peaks are identified around some frequencies. Based on the 55 
days of photometric observations alone, it is hard to make the 
distinction between real multiplets and non-stationarity of the 
amplitudes. 

Close to the 8 Set instability strip, there is another 
small group of multiperiodic variables: CoRoT 102933855, 
102703484, 102850502, and 102816758. Although these stars 
have a rich frequency spectrum, they are different from the pre- 
vious group because their spectrum shows more isolated fre- 
quencies with higher amplitudes (~ 0.5%). Schematic oscilla- 
tion spectra are provided in Fig. |A.5l 

Near the border of the SPB instability strip, four monope- 
riodic variables are visible: CoRoT 102921797, 102774512, 
102889144, and 102872474. The distinction between monope- 
riodic oscillation and spots is difficult to make without colour or 
spectral information. 
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Fig. 8. Selection of spacing test results: (left) the inverse variance test of detected frequencies in CoRoT 102922479 indicates a 
spacing of 0.78 d . The matrix number (middle) test gives a clear spacing of 0.74 d -1 in CoRoT 102848985. The Kolmogorov- 
Smirnov (right) test applied to CoRoT 102833548 reveals a spacing of APi = 0.66 d , but does not exclude APi /2 as an alternative 
possibility. 
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Fig. 9. Time-frequency analysis of CoRoT 102822479 (left), 102848985 (middle), and 102833548 (right). The spacings found (see 
Fig. |8]l and their derived frequency members are depicted by the horizontal short lines, vertically connected. 



The main power of all of the above stars is located in the ex- 
pected (3 Cep pulsation range. However, some of the stars on the 
red side of the SPB instability strip show variability on longer 
time scales (< 4d~'), e.g. Figs lA.6l and IA.7I in Appendix A. 
Not all of this periodic variability is sinusoidal or even approx- 
imately sinusoidal. The best example in this respect is CoRoT 
102762284, for which the sinusoidal variation is superposed on 
a larger scale outburst-like variation (Fig. IA.81 l. 

A summary of the interpretation of the cool, short-term (i.e. 
P Cep-like) candidate pulsators in the log 77 e ff - log g diagram is 
shown in Fig. [12] 

4.3. SPB candidate examples 

As is clear from Fig. [71 the SPB star candidates from the initial 
classification agree with their instability strip. The 15 clearest 
SPB candidates, in order of decreasing amplitude, from ~ 1 % 
to ~ 0.2%, are CoRoT 102844894, 102754851, 102728830, 
102855391, 102739246, 102797587, 102938439, 102769848, 
102848506, 102848506, 102956197, 102943966, 102968674 
and 102871668. Their effective temperature and gravity locates 
them well within the theoretical SPB instability strip. Besides the 
typical closely spaced g mode pulsations, we also find a resid- 




encies are phase-binned data points, the light grey line is a fit us- 
ing 7 harmonics, determined via consecutive prewhitening. 



ual power excess bump in a frequency band above 1.5-2 d for 
the majority of these stars, well separated from the typical SPB 
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Fig. 10. B-type pulsator candidates on the red side of the SPB instability strip. Some of these stars show clear, stable oscillations. 
(top) CoRoT 102816758 (f =2593.0705 CJD, /, « 3.6d _1 ), (middle) CoRoT 102833548 (f =2593.0705 CJD, f x « 4.0d _1 ), (bot- 
tom) CoRoT 102850502 (fo=2593.0705 CJD, f\ » 5.0d~'). Light curves with high time resolution were binned per ~10 datapoints 
for visibility reasons. The black line represents the fit from which the satellite orbital frequency, its harmonics, and its aliases were 
filtered out. 
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Fig. 14. (top panel) Short Time Fourier Transformation of the residual light curve of CoRoT 102938439, after prewhitening a 
nonlinear fit with all detected significant frequencies below ~1.5d _1 . The highest amplitude frequencies in the residuals clearly 
show modulated amplitudes, (bottom panel) A binned version of the residual lightcurve clearly shows the largest temporal variability 
near day 25 of the light curve, with some smaller 'outbursts' as well. 
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Fig. 12. Location of (3 Cep candidate pulsators in the log T e g — 
log g diagram with a determined T e g too low to be a real f3 Cep 
star: 5 Set like variables (black dots), monoperiodic (3 Cep can- 
didates (white dots), rich pulsators (grey diamonds) ,and targets 
with equidistant frequency spacing (black white diamonds). 
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Fig. 13. CoRoT 102938439: excerpt of original light curve and 
fit (solid grey line), around a time where the amplitude is tem- 
porarily driven to a higher level. The typical shape of shallow 
minima and sharp maxima is visible. 

pulsation range predicted for excited modes of £ < 3 modes. 
The theoretical instability strip reaches the observed range of 
frequencies for modes with a higher t value, but at the same 
time, disk-averaging effects reduce the observed amplitudes. For 
t = 4, this reduction factor amounts to about 30, which is still 
above the CoRoT noise level for many of the SPB candidates, 
so we cannot exclude the higher observed frequencies being due 
to high-degree modes. We illustrate the nature of the observed 
variability for a prototype of these stars (CoRoT 102938439) in 
Figs.[14j |T3] and Qj] Other examples can be found in Figs. |A.9l 
to |A.12t . After prewhitening all low frequencies, the largest 
feature in the amplitude spectrum turns out to be a collection 
of closely spaced peaks (Fig. [15]). At least two phenomena can 
cause these features: unresolved stable modes with a complex 
beating pattern, or time-dependent oscillation frequencies and/or 
amplitudes. In the first case, we expect a region with blended or 
nearly blended frequency peaks, without any amplitude struc- 
ture. In the second case, we expect a frequency region with 
a peak at the excited frequency, and Lorentzian shape farther 



0.6 




Fig. 15. Original periodogram of CoRoT 102938439. The inset 
is the amplitude spectrum of the residuals, after prewhitening 
all significant frequencies below ~ 1.5 d _I . Grey solid line is a 
Gaussian smoothing with cr = 0.075 d -1 . The black solid line 
denotes the mean noise level, the dotted line is the 99% signifi- 
cance level. 



away from the central peak. When a residual light curve of these 
stars is constructed by prewhitening a nonlinear least squares fit 
of all frequencies below ~1.5d _1 , the same amplitude structure 
emerges for all of them, which points towards the same physical 
origin given that the pattern is not chaotic. The amplitudes of 
the detected modes are highly variable over time (see Figs |A.9l 
and lA.lOl i. The observed features could come from modes with 
a finite lifetime or to nonlinear nonresonant distortion leading to 
time-dependent phenomena. An analysis of the evolution of the 
amplitudes shows that the period over which the modes appear 
and dissappear is of the order of one to a few days. Instead of a si- 
nusoidal shape, the light curve maxima are narrower and sharper, 
while the minima are shallow and broad (Fig.[l~3l. This is very 
similar to the light curves of some white dwarf pul sators and 
is exp lained by nonlinear effects for these objects (Vuill e~et"aT1 
2000). We notice that a similar feature of disappearing and reap- 
pearing frequencies was also found in the Be star HP 49330 an d 
was interpreted as the cause of an outburst dHuat et al.l 12009). 
It could very well be that this is a general phenomenon at low 
amplitudes for B stars, but that it only leads to mass loss for the 
most rapid rotators among them. 

A second manifestation of nonlinearity can be found in 
many detected combination frequencies. In contrast to the time- 
dependent phenomenon described above, these are stable on the 
scale of the total time span, thus resulting in isolated frequency 
peaks in the spectra (see online frequency tables). A detailed 
description of these modes is difficult based on the short time 
span of the Initial Run and the long periods of these modes 
(< 2 d~'), but will surely be possible for similar new pulsators in 
the CoRoT long run data. 

Secondary good SPB star candidates include CoRoT 
102863407, 102752912, 102773435, 102917802, 102826973, 
102838201, 102888003, 102915048, 102813396, and 
102945383 of which the last four also show evidence for 
modes of finite lifetime. 



10 



P. Degroote et al.: Asteroseismology of B stars in CoRoT's Initial Run Exofield 



4.4. Other variable objects 

Brightness variations need not necessarily be due to pulsa- 
tions (alone). For example, when two binaries are physically 
close to each other and the stars are of comparable luminos- 
ity and size, and the grazing eclipses can resemble one single, 
stable sinusoidal pulsation. A DFT analysis puts more power 
in the double frequency of the orbital frequency, folding both 
the primary and secondary eclipse onto each other. Therefore, 
close binaries with an orbital period roughly between 0.25 and 
0.5 d~ l can easily be mistaken for a ySCep star, whereas a 
wider system can be misinterpreted as an SPB. Examples of 
close or contact binaries are CoRoT 102725806, 102858055, 
102876625, 102930503, 102784048, 102768686, 102732139, 
102873761, and 102812372. The phase diagram of the first can- 
didate, CoRoT 102725806, is shown in Fig.QT] We find an or- 
bital frequency of / = 2.925(7) d~', independently confirmed 
with the phase dispersion minimization and information entropy 
method. The amplitude of the primary eclipse is A\ = 8.56%, 
compared to the average flux level. The double system was mis- 
taken for a single /?Cep star by the CVC because of the small 
difference between the amplitude of the primary and secondary 
eclipse; A\-Ai = 0.87%. The duration of the primary eclipse in 
phase units is Pi « 0.08, which is almost equal to the duration 
of the secondary eclipse, P2 ~ 0.10. 

The possibility of a binary system and an SPB candidate are 
not mutually exclusive, because CoRoT 102918586, 102912741 
and CoRoT 102793963 turn out to be candidate SPBs in close 
binaries, of which the latter shows a single sinusoidal variation, 
synchronised with the orbit of the binary. These systems in par- 
ticular can be of interest for asteroseismology, because the bi- 
nary orbit can help to constrain the mass and radius of the SPB 
star. 

Another group of SPB candidates have frequencies in a very 
broad range, e.g. CoRoT 10282445, 102764403, 102818535, 
102816987, 102837646, 102804931, 102887852, 102930369, 
102790135, and 102856178. CoRoT 10282445, for example, 
shows periodic variations with time scales around - Id -1 , 
~4d 1 and ~15d~' (Fig. [TBI. The main frequency f\ 
1 .04 1 1 (7 ) d can be connected to the third frequency fa = 
1.1533(3)d~', resulting in a beating period of almost 10 days. 
The second frequency fi = 15.271(9) is connected to fa = 
14.739(4)d~', and the same spacing of 0.52d _1 is found 
throughout the whole frequency spectrum (Table lATl ). Although 
the spacing connects the three main frequency regions where 
variability occurs, current stellar models do not offer any inter- 
pretation in terms of a single SPB, /3 Cep, or 6 Set star, while we 
can exclude contamination from nearby stars on the CCD. The 
possibility of a binary system consisting of a hybrid /? Cep and a 
6 Set remains, as well as the one of a triple system. Spectra will 
help us to unravel the exact cause of the variation in this light 
curve. 

For a couple of light curves (CoRoT 102910610, 102968526, 
102857565, and 102876625, see Fig.[l7]l, spot models seem 
more appropriate for describing them than a pulsation model; the 
resemblance to simulate d light curves of rotating spotted stars 
is ind eed apparent (e.g. Strassmeier & Boppl [19921: lLanza etakl 
1 1993b . Their frequency spectra agree with those models, and 
few frequencies are detected in a very narrow interval. Figure[T8l 
illustrates that it can be difficult to distinguish simple beating 
patterns from spot related variability. By assuming that the star 
rotates differentially like the Sun, the multiplet structure can be 
used to determine the level of differential rotation AO., e.g., for 
CoRoT 102968526, AQ = 0.0538 +0.0002d-\ This implies that 




Fig. 19. Theoretical instability domains of £ — 0, 1, 2 modes rep- 
resented in a log r e (j-log P diagram for main-sequence models 
with Z = 0.02 (upper panel) and Z = 0.01 (lower panel). In each 
panel, the regions of unstable modes represent /3 Cep (a), SPB 
(b), S Set (c), and y Dor-type pulsations (d).The frequencies de- 
tected in the best SPB candidates (see text) and the candidate 
pulsators from Sect. l4.2l are shown as triangles. 



one part of the star completes one more full rotation cycle over 
a period of 1 / AQ. 

5. Comparison with theoretical predictions 

The results presented in the previous section show that many of 
the stars classified as an SPB or /3 Cep star by the CVC classi- 
fier are not typical members of the two classes (at least as far 
as the classes have been defined so far). To have a general com- 
parison between observed and theoretically predicted instability 
domains, we not only consider log T e ff - \ogg planes (see Fig. 
13, but also confront predicted and observed frequency instabil- 
ity domains, as presented in Fig. [19] 

The theoretical instability domains were computed for non- 
rotating main-sequence models of masses between 1.2 and 18 
M and 2 values of the heavy-elements mass fraction: Z=0.02 
and Z=0.01 (respectively upper and lo wer panels of Fig . [T9l . 
Stellar models were computed wit h cles (IScuflaire et al.l2 008b). 
adiabatic frequencies with losc dScuflaire et alj|2008al) . and the 
stability of mo des of degree t — 0, 1,2 was investigated with 
the code mad dDupret et al.l 120031) . In the case of S Set and 
y Doradus stars, the nonadiabatic computations included the 
interaction between con vection and pulsation as described by 
Grigahce ne et al.ld2005l) . All stellar models were computed with 
overshooting from the convective core (aov = 0.2) and a 
mixing-length parameter <?mlt = 2.0, both expre ssed in local- 
pressu re, s cale-height uni t s. OP opacity tables dBadnell et al.l 
120051) with lAsplund et all d2005l) met al mixture were adopted. 
We re fer to Miglio et al. (200 7albl) and lZdravkov & Pamvat nvkh 
(2008) for a detailed study of the effect of considering different 
opacity tables and metal mixtures on the excitation of pulsation 
modes in B-type stars. The distinct islands of instability populat- 
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Fig. 16. Phase diagrams of three clear sinusoidal variations of CoRoT 102804522. The frequency in the left panel (/ a 1.04d _1 ) 
is rather high for typical g modes in SPBs, while the frequency in the middle panel (/ ss 15.27 d _1 ) is rather high for a typical 
ft Cep star p-mode, but quite typical for p modes in 6 Set stars. The frequency with which the light curve is folded in the right panel 
(/ « 4.05 d _1 ) could, on the other hand, be a typical p mode in a /3Cep star. This different behaviour raises the suspicion that this 
is not one single star, although there are neither traces of eclipses nor any indication of flux contamination from nearby stars on the 
CCD. 




17 

t - i (CoRoT JD) 

Fig. 17. Selection of ySCep candidates showing strong hints of spot features in their light curves, (upper panel) CoRoT 102968526 
(?o=2590.0450 CJD) has a dominant frequency of /i * 4. 14 d -1 with a secondary peak at fa * 4.08 d _1 . The black line is a fit using 
all noninstrumental frequencies (bottom panel) The light curve of CoRoT 102857565 (fo=2590.0499 CJD) was binned per 10 points 
for visibility reasons. 



ing the log T^s -Period diagram of Figs[T9l correspond (from left 
to right) tojSCep-, SPB-, 6 Set-, and y Dor-type pulsations. 

Stars classified as SPBs are in overall agreement with the 
theoretically expected domain (see Figs. [7] and 119) , though in 
some cases with frequencies higher than those of a typical SPB 
star (which could be due to the strong effects of rotation on the 
periods of oscillation). On the other hand, several stars classified 
aSjSCep are in the typical T e g instability domain of SPBs. These 
stars have frequency peaks in the f3 Cep frequency domain: their 
modes could be p rograde SPB-type modes as proposed in, e.g., 
ISaioetalJ d2007l) . 

The effects of rotation may also affect the determination of 
T e ff for a few stars that have a frequency spectrum resembling 
the one of fast-rotating SPBs but that, based on their position 



in Fig.|T9j are cooler than the red edge of the SPB instability 
strip (e.g. CoRoT 102833548 and 102729531). A more reliable 
spectroscopic r e ff determination of these stars is thus needed to 
clarify the nature of these targets. Such spectroscopic data will 
become available in the next months. 

The comparison with theoretical expectations will also 
greatly benefit from the determination of the number of pulsat- 
ing B stars relative to all the B stars observed in CoRoT's IRaOl 
EXOfield. The absence of typical /3Cep pulsators (in terms at 
the same time of frequency spectrum and T e ff) can be ascribed 
simply to the expected scarcity of early B-type stars in the field, 
following the initial mass function and stellar evolution models. 
This requires the determination of stellar parameters for a large 
number of stars and will be addressed in a future work. 
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Fig. 18. CoRoT 1028355 15 (upper panel) and 102824452 (lower panel, light curve binned per 10 points) have a frequency spectrum 
where the power is concentrated in narrow bands at multiples of the dominant frequency. Combined with the peculiar shape of their 
light curves, they strongly resemble spotted stars, although the distinction from a simple beating pattern is not easy (bottom). The 
solid line is a fit, grey circles are data. 



It is clear from these first general remarks that the relevance 
of CoRoT's photometric variability survey goes beyond the sim- 
ple definition of instability domains of expected and known 
types of pulsators. Particular theoretical interest is provided by 
the unexpected variability detected in several targets on the red 
side of the SPB instability strip, a region where several claims 
of the detectio n of pulsations were made (the so-called "M aia" 
stars, see e.g. lScholz et al ] |!998tlAerts & Kolen berg 2005|). The 
theoretical instability mechanisms proposed as an explanation 
for these stars do not cover the full area in the Hertzsprung- 
Russell diagram wh ere such pulsati ons have been claimed (see 
lTownsendll2005l and Savoniie 2005). The forthcoming spectro- 
scopic observations of those targets will shed new light on these 
most interesting targets, and allow theoretical interpretation of 
the detected variability. 

6. Conclusions 

We constructed a frequency analysis pipeline and additional 
post-processing tools to efficiently analyse pulsators in CoRoT's 
exoplanet database. The method also includes a jump correction 
and detrending algorithm, as well as an automated search for 
combination frequencies and spacings among frequencies or pe- 
riods. In the particular case of the Initial Run, many combination 
frequencies and spacings were found, but the total time span of 
~ 55 d was too short to draw firm conclusions on the frequency 
or period spacings of g modes and to confront them with the- 
oretical models. There were too few convincing candidates in 
the more massive /3 Cep class to commence an interpretation in 
terms of nonlinear resonant coupling. However, our techniques 
revealed clear frequency spacings in high-frequency variables on 
the cool side of the SPB instability strip. 

Most of the stars that seemed good /3 Cep candidates in 
terms of the values of the first three frequencies are to be re- 
futed. This was to be expected, because the number of can- 
didate /?Cep stars foun d by the CoRoT Variability Classifier 
dDebo sscher et al. 2009) was anomalously high compared to the 



number of candidate SPBs. Only one star on the blue side of the 
SPB instability strip matches a typical /?Cep frequency spec- 
trum. Four other stars were found to be 6 Sct/yS Cep-like, but 
again were far off their theoretical location in the log r e (j - log g 
diagram. Instead, we detected many different types of variable 
stars among the ySCep candidate sample. Several turned out to 
be low-amplitude pulsators whose frequencies are spread over 
a wide range and which are located on the red side of the SPB 
instability strip, where no oscillations are pr edicted by t heory - 
Acco rding to the th eoretical computations bv lTownsendl ([2005) 
and lSavoniid (120051) . mid-B type rapidly rotating pulsators with 
retrograde mixed g modes could occur in a part of this region, 
exhibiting many of these properties. However, they could also 
simply be fast-rotating classical SPB or 6 Set stars, or Be pul- 
sators. Future spectroscopic observations will allow us to distin- 
guish between these possibilities, by supplying information on 
rotational velocities, abundances, and improved log g - log 
values. Besides these interesting stars, also binary stars and a 
few suspectedly differentially rotating stars (with or without ad- 
ditional pulsations), were detected among the /? Cep candidates. 

Similar conclusions apply to the sample of SPB candidates, 
although the number of good candidates is vastly greater for this 
class than for the ySCep candidates. Interestingly, the highest 
amplitude SPB candidates all show two distinct signs of non- 
linear behaviour: first, clear unexpected power excess in the re- 
gion between 2 and 6d~\ which seems to point to modes with 
finite lifetimes. These could originate from nonlinear nonreso- 
nant harmonic distortion, leading to time-dependent amplitudes 
and/or frequencies. Second, we identified combination frequen- 
cies, which are stable on the time scale of the measurements. 
These are likely to be caused by nonlinear resonant mode lock- 
ing. On top of that, we also find some rich and low-amplitude 
SPB pulsators, with frequencies in a lower range than the simi- 
lar stars in the /? Cep sample. 

Our follow-up studies of new pulsating B stars in the CoRoT 
long runs, which typically last 150 d, as well as future spec- 
troscopy, will be very valuable in evaluating the new type of os- 
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cillations found in the Initial Run data as reported here in more 
detail. 
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Fig.A.l. Sames as Fig. IA.4l but for the stars showing a com- 
plex but structured frequency spectrum (CoRoT 102833548, 
102848985, 102922479, 102850576). 
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Fig. A.2. /3 Cep candidate 102729531: {bottom panel) Zoom on day 17 of IRaOl: on top of the clear ~ 0.25 d period, many more 
significant variations of the flux are noticable. {upper two panels) raw frequency spectrum (grey) with identified frequencies with 
SNR> 4 (black lines) and orbital frequencies (grey lines, all harmonics and day aliases of 1 to 5 days). 
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Fig. A.3. Same as Fip.DOl but for the /3 Cep candidate 102790063. 
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Fig. A.4. Schematic oscillation spectra of one group of pulsators 
of which the fundamental parameters log T e ff and log g place 
them in between the SPB and 5 Set instability strip (CoRoT 
102861067, 102790063, 102729531, 102862454, 102790331, 
102771057). These are all multiperiodic variables, showing 
complex variability in a broad frequency range (independent 
modes are depicted with black solid lines, harmonics as inter- 
rupted black lines, and candidate combination frequencies in 
grey). 
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Fig.A.5. Same as Fig. IA.4l but for two of the stars with a rich 
frequency spectrum, but having few large-amplitude pulsations 
(CoRoT 102850502, 102816758). 



0.0010 



0.0005 



< 



-3 o.oooo 



-0.0005 




0.4 0.6 
Phase 



Fig. A.7. CoRoT 102808565: Binned phase diagram (lcr errors 
on average), after removal of a nonlinear fit with 69 frequencies. 
Inset: excerpt from the residual Scargle periodogram. 
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Fig.A.8. CoRoT 102762284: Phase diagram of the residual 
lightcurve after removal of the first two frequencies, which are 
both well approximated by sinusoidal variations. The diagram 
was folded with frequency / = 0.278 12 d~'. A smaller period 
variation is noticable, as well as an outburst-like peak in the sec- 
ond half of the phase. 
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Fig.A.6. CoRoT 102872265 {upper panel) and 102832493 (lower panel) show a complex frequency spectrum with hig 
with frequencies above the typical SPB pulsation range. Moreover, the dense spectrum makes it difficult to seperate lonj 
pulsations from beating patterns. The solid line is a fit, grey circles are the data. 
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Table A.l. CoRoT 102804522 shows rich variability on vary- 
ing timescales, but many of the frequencies can be connected 
through the same frequency spacing (e = fa — f\- A/). 



Type 


h (d- 1 ) 


A/(d-') 


/2(d-') 


<f(cT) 


Single 


1.0891 


0.5320 


1.6213 


0.0002 


Single 


1.2350 


0.5248 


1.7595 


-0.0003 


Single 


1.4445 


0.5248 


1.9697 


0.0004 


Single 


1.7595 


0.5320 


2.2914 


-0.0002 


Single 


1.8758 


0.5320 


2.4077 


-0.0002 


Single 


3.1081 


0.5248 


3.6325 


-0.0004 


Single 


3.6325 


0.5248 


4.1572 


-0.0001 


Single 


3.7455 


0.5248 


4.2703 


-0.0001 


Single 


3.8578 


0.5320 


4.3899 


0.0001 


Single 


3.9212 


0.5248 


4.4458 


-0.0003 


Single 


4.3259 


0.5248 


4.8512 


0.0004 


Single 


4.7765 


0.5320 


5.3088 


0.0003 


Single 


9.0494 


0.5248 


9.5741 


-0.0001 


Triplet 


14.2073 


0.5320 


14.7395 


0.0001 




14.7395 


0.5320 


15.2719 


0.0004 
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Fig. A.9. CoRoT 102943966: Original periodograms of SPB candidates showing residual power excess at frequencies above ~1.5 
d _1 (upper left). The insets are amplitude spectra of the residuals, after prewhitening all frequency in lower regions. The grey solid 
line is a smoothed version of the spectrum, horizontal lines denote theoretical noise level and 99% confidence interval, upper right: 
Short time Fourier transformations of the residual light curves give the evolution of the (normalised) power in the periodogram over 
time since the start of the measurements (t). This time-frequency analysis shows that the excess power is due to modes with a finite 
lifetime in the order of days, (lower panel:) excerpt from the light curve, showing the shape of the flux changes due to the pulsations. 
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Fig. A.10. CoRoT 102855391: same as Fig.lA^ 
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Fig. A.ll. {left column) Original periodograms of SPB candidates showing residual power excess at frequencies above ~1.5 d _1 . 
The insets are amplitude spectra of the residuals, after prewhitening all frequency in lower regions. The grey solid line is a smoothed 
version of the spectrum, horizontal lines denote theoretical noise level and 99% confidence interval, (right column) excerpt from the 
lightcurve, showing the shape of the flux changes due to the pulsations, (from top to bottom) CoRoT 102769848, 102754851 and 
102848506. 
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Fig. A.12. Same as FigjAJj] but for CoRoT 102844894 and 102956197. 



